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1 Introduction 



The mass of the tau-neutrino is very loosely bound by experiments [Q: 

m^^<18MeV, (1) 

whereas it is much stronger restricted by cosmology. In the case that the Vj. is stable 
on the scale of the universe age, tu ~ 10 Gyr, its mass obeys the Gerstein-Zel'dovich 
limit [0, i.e. roughly < 10 eV. If it is unstable, the cosmological limit is much 
less restrictive and depends upon the life-time of Vr and possible decay channels. 

If the Vr is unstable on the cosmological time scale, so that the Gerstein-Zel'dovich 
limit is avoided, but is stable on the nucleosynthesis time scale, i.e. r^,^ > 200 sec, 
then considerations of primordial nucleosynthesis lead to a much better bound than 
(|I|). The bound on the i/^-mass from nucleosynthesis was first found in ref. 0] and 
in a slightly improved form in ref. [Q. The calculations of the second work predict 
a somewhat larger value of the frozen energy density of z/^-, but in the translation of 
this result to the effective number of neutrino species, found from the distortion of 
^ife abundance, a numerical error was done which resulted in an overestimation of 
the number of additional effective neutrino species. Still, even with the correction of 
this error the results of ref. ||^ are stronger than those of the pioneering paper [^. 
The calculations of both papers were done under the following basic assumptions. 
It was assumed that the massive and the two massless neutrinos, and i/^, are 
in complete kinetic equilibrium so that their energy distributions are given by the 
canonical expressions: 

exp{(3E-Q + 1 

where /3 = 1/T is the inverse temperature and are dimensionless effective chemical 
potentials. Another two simplifying assumptions were made, namely that the chemi- 
cal potentials of massless neutrinos are zero, C,^^ = (,u^ = 0, and that the distribution 
function of Ur can be approximated by its Boltzmann limit: 

/(E) = exp[(^(t)-i5;)/T(t)]«l, (3) 



which is accurate when the temperature is small in comparison with the mass, m > T. 
Here we use the standard notation for the chemical potential, /i = ^T. 

It is well known that with these approximations the system of complicated integro- 
differential kinetic equations (see eq. @ below) is reduced to a single ordinary differ- 
ential equation (eq. (^) for the number density of Uj. or, equivalently, for the unknown 
function ^(t) (see e.g. books ^j), which easily can be solved numerically. A similar 
approach was recently used in ref. where the role of a massive Ur in the production 
of all primordial light elements and not only of '^He, as in refs. ^, was considered. 

The assumption of vanishing C,ue was relaxed in ref. If in addition 

one assumes validity of Boltzmann statistics, then instead of one differential equation 
for there appear three ordinary coupled equations for three chemical potentials 
C,j{t),j = e, fi, T. If one uses exact Fermi-Dirac statistics, then the unknown functions 
^j(t) remain "inside" the collision integrals but the kinetic equations still remain 
ordinary differential equations in this approximation, and not integro-differential ones 
for fj{t,p) as in the exact case. This approach permits one to take into account an 
average heating of massless neutrinos by annihilation of massive tau-neutrinos. It 
naturally results in a smaller neutron-to-proton ratio and in a weaker influence of a 
possible non-zero m^^ on primordial nucleosynthesis. On the other hand, one can 
easily see that nonzero C,ue and ^j^^ (in kinetic equilibrium) give rise to a considerably 
larger frozen number density of z/^- Correspondingly the net effect on nucleosynthesis 
is stronger. We do indeed obtain a large increase in n^^, however, a x^"fit, for different 
Ur masses, of our resulting distribution functions for the massless neutrinos to exact 
FD-distributions allowing for both an effective temperature (possibly different from 
the expected neutrino temperature) and a chemical potential gives C,^^ in the range 
from -|-10~^ to —2 ■ 10^^ with extremum near m^^ = 5MeV. 

It was found in ref. [§] that together with an overall neutrino heating there are 
considerable distortions at the high energy tail of massless neutrino spectra, created 
by annihilation of heavy tau-neutrinos into massless (or light) z/'s. The spectral 
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distortion is especially important for the electronic neutrinos because they directly 
influence the frozen n/p-ratio, through the reactions: 

UeU ^ pe~ and e^n <-> pue . (4) 

The calculations of the spectral distortions were made in the Boltzmann limit and the 
problem was reduced to the solution of an ordinary differential equation. According 
to ref. IP the deviation from kinetic equilibrium of the electronic neutrinos is quite 
essential and has a noticeable influence on nucleosynthesis. These calculations were 
refined in [0. The calculations of the present paper are in agreement with the semi- 



analytical estimates of the distortion of Ue spectrum made in refs. |10[, which results 
in a larger n/p-ratio. The spectral distortion is not just a peak from z/T--annihilation, 
since the excessive energy is redistributed by elastic scattering over the whole spectra 
of z/e,^. This effect gives rise to an overall neutrino heating or, in other words, to 
larger temperatures of the massless neutrinos relative to the photon temperature. 
This leads in turn to a later freezing of the reactions (^) and to a lower n/p-ratio. 
According to our calculations the average heating produces a stronger effect than the 
spectral distortion. This agrees with the statement made in ref. (see also "Note 
added" to ref. §). 

In view of this discussion it is interesting to solve the exact system of integro- 
differential kinetic equations numerically (with a good precision) and to find exactly, 
without any simplifying assumptions, the impact of possibly massive tau-neutrinos on 
primordial nucleosynthesis. The first calculations of this kind were done in ref. [|lT 



In what follows we will correct some matrix elements for scattering of Majorana 



neutrinos presented in ref. [|lT|, and solve the kinetic equations numerically, using a 



more accurate code. Our results for n,^^ are in a good agreement with ref. |jTT[. We also 



agree with this paper in the calculation of the changes in light element abundances for 
low values of rriu^. In the high mass range our results are noticeably lower. We ascribe 
this difference to a lower momentum cut-off made in ref. [|TT], i/max ~ ip/T)max ~ 13, 



in the region being essential for the reactions (§), while we made the calculations with 



a considerably larger cut-off, i/max = 20. 

Our technique repeats that of ref. |jl2[, where similar calculations were done for 



massless neutrinos, with evident complications in the case of a massive due to 
a larger number of equations and different matrix elements. We conclude that the 
influence of non-equilibrium corrections on primordial nucleosynthesis is considerably 
larger than what was calculated in the earlier papers under assumptions of kinetic 
equilibrium and Boltzmann statistics. The dominant effect comes from an increase 
of the frozen energy density of massive z/,-. 

2 Qualitative discussion 

The complete set of integro-differential kinetic equations has the form: 

{dt-Hpd,)f){p„t) = lf\ (5) 

where the collision integral for two-body reactions 1 + 2 ^3 + 4 is given by the 
expression: 

rcoU ^ J_ V f d^P^ d,^P3 d:^P^ n\M2 

^ 2Ei ^7 2E2(27r)3 2E3(27r)3 2E4(27r)3 ' '12-^4 

(27r)^5(^)(pi +P2-P3- p,)F{f„ /2, /3, U) , (6) 

where F = fsf^il - /i)(l - - /i/2(l - /3)(1 - Ja), I^P is the weak interaction 
amplitude squared summed over spins of all particles, and S is the symmetrization 
factor which includes 1/2 from the averaging over the first particle, 1/2! for each 
pair of identical particles in the initial and final states and the factor 2 if there are 2 
identical particles in the initial state; the summation is done over all possible sets of 
leptons 2, 3, and 4. 

Numerical solution of such a set of equations is tremendously more difficult than 
the solution of the well known simple ordinary differential equation for the number 
density of massive neutrinos, which is valid in the case when the neutrinos are in 
kinetic equilibrium and obey Boltzmann statistics: 



riur = {crannv)in^g-n^^). (7) 



One would naturally ask if such a complexity is necessary. Eq. gives a very good 
approximation if the rate of elastic scattering of the massive particles in question is 
much larger than their annihilation rate. This is usually true for very non-relativistic 
particles since their annihilation is Boltzmann suppressed while elastic scattering 
is not. For the case of z/^, both elastic scattering and annihilation are frozen at 
temperatures which are not much below m^^ and rather close to each other so the 
effects of deviations from equilibrium may be significant. When one calculates the 
frozen number density of massive Ur, a gradual switch-off of elastic scattering results 
in a distortion of their spectrum and in a faster-than-equilibrium cooling. Since the 
cross-section of their annihilation is proportional to the energy squared, this faster- 
than-equilibrium cooling gives rise to an increase in their frozen number and energy 
densities and results in a larger effect on nucleosynthesis. A simple attempt to take 
this non-equilibrium cooling into account was made in ref. ^ where it was assumed 
that at temperatures above a certain value, Teq, which is determined by the strength 
of the elastic scattering, tau-neutrinos have the equilibrium distribution in energy. 



f,y^ = exp[— + ml^/T + ^], while below T = T^q the complete kinetic decoupling 
was assumed, so that the momentum p was scaled as inverse universe expansion 
factor a{t). Unfortunately this simple anzats is not very accurate and more elaborate 
calculations (ref. ||l3l) show, that the spectrum is distorted in a more complicated way. 



Rough estimates of paper [jT3[ show that these non-equihbrium corrections give rise 
to an increase of n^^ at the level of 10%. As we see below in fig. 1, the calculations of 
the present paper show an even bigger effect. The frozen number density of heavy 
differ from that of ref. by 35 % and from that of ref. by 15 % in the maximum. 

Another important effect is a distortion of the equilibrium spectra of massless 
neutrinos and in particular of electronic neutrinos. An excess of z/g at the high energy 
tail of the spectrum gives rise to a larger frozen number density of neutrons and corre- 
spondingly to a larger abundance of '^He. There is also an opposite effect related to an 
overall increase in the number density of electronic neutrinos due to the annihilation 
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u^ii^ — i> z/eZ/g. This results in a later neutron decoupling and to a smaller frozen ratio 
n/p = exp(— Am/Tj), where Tf is the freezing temperature of the (n ^ p)-reactions 
and Am = 1.293 MeV is the neutron-proton mass difference. It is difficult to sep- 
arate rigorously these two effects using the z/g spectrum calculated here, because of 
an ambiguity in the choice of the unperturbed distribution, in particular because the 
effective temperature is different in different parts of the spectrum. 

It is also important that the ratios of the photon temperature and the (average) 
temperatures of the massless neutrinos are different from the canonical value T^/T^ = 
(11/4)^/'^. It can easily be understood that these ratios depend non-monotonically on 
m^^. They go up with m^^, reach maximum near m,y^ = 5 MeV and fall down close 
to the massless value for large rriy^. 

One sees that there are several essential effects associated with the distortion of 
kinetic equilibrium which significantly may change the primordial abundances of ^ife 
and other light elements. All of them are automatically taken into account in our 
calculations below. 

3 Kinetic equations 

We want to solve kinetic equations numerically for the following system of interacting 
particles: massless z/g and u^, massive z/^, electron-positron pairs, and quanta of 
electro-magnetic radiation. There are three unknown distribution function of time 
and momenta, fj{t,p), for the three neutrino species. We assume that electron- 
positrons and photons are in complete equilibrium so that their distributions are given 
by eq. with vanishing chemical potentials and with temperature T{t), which is an 
unknown function of time (of course photons obey Bose statistics). The assumption of 
equilibrium of the electro-magnetic component of the primeval plasma is very accurate 
because of a large strength of the electro-magnetic interactions. Thus there are 3 
unknown functions of t and p and one unknown function of t. For their determination 
we have three kinetic equations (|^). The fourth necessary equation is the covariant 
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conservation of energy- momentum in the expanding universe: 

p = -SHip + P) . (8) 

Usually one determines the temporal evolution of temperature from the law of entropy 
conservation. However, the latter is only true in thermal equilibrium, so if one wants 
to describe deviations from equilibrium consistently, a more general, though more 
complicated eq. (H) must be used. In this equation p and P are respectively total 
energy and pressure densities in the cosmic plasma. They are given by the expressions: 



P- 15 +ivr2[exp(E/T,) + l]+^.i:,^i vr^ +i ,2 + (9) 



and: 



45 ~^J3n^[exp{E/T,) + l]^^£^J 3rr^ J Srr^J^^^^ ^ ' 



The Hubble parameter, H = d/a, is related to the total energy density p in the usual 
way: 



inp 
^pi 

119 



= , (11) 
Sm, 



where a{t) is the universe expansion factor (scale factor) and mpi = 1.22 ■ 10 GeV 
is the Planck mass. 

It is convenient instead of time and momenta to use the following dimensionless 
variables: 

X = ma{t), Uj = Pja{t), (12) 

where m is an arbitrary parameter with dimension of mass, which we took as m = 
1 MeV, and the scale factor a{t) is normalized so that a{t) = l/T^, = l/T^ at high 
temperatures or at early times. In terms of these variables the kinetic equations (^) 
can be rewritten as: 

HxdJ,{x,yr)=lfK (13) 

The relevant reactions and the corresponding matrix elements squared are presented 
in tables 1 and 2 for the cases when the first particle is iy^) and respectively. 



There are the following differences with ref. |]TT|: 1) for the reactions Uai^a i^ai^a we 
took twice larger contribution because of identical particles in the initial state; 2) we 
have a different expression for the matrix element squared for elastic scattering of 
massive Major ana Vr- 

Our procedure of solving these equations is essentially the same as in our previous 
paper [T^, where non-equilibrium corrections to the spectra were calculated for the 



case when all neutrinos were massless. The collision integral is reduced from 9 down to 



2 dimensions by the method described in ref. |T^ with some complications connected 
with the momentum dependence of the matrix elements (see ref. |12|)- The results of 
the calculations are presented in Appendix A. 

3.1 The problem of initial conditions 

We solve the system of kinetic equations (pTsD in the "time" interval Xj„ < x < Xout- 
All the collision integrals in the r.h.s. of eq. ( |T3| ) are suppressed at large x at least 
by the factor We find that at x ~ 50 all the variables, we are interested in, 

reach their asymptotical values so, to be on the safe side, we choose the final time 
Xout = 100. 

The nucleosynthesis code of ref. [|15| requires the final time x ~ 2000 for all energy 



densities and n ^ p rates. We applied a separate program, which calculates these 
quantities, using the energy conservation law and our final values for the distribution 
functions /y- and the temperature T^. 

At early times, x ^ 1, the collision rates in eq. (|TB|) are very high due to the 
factor and therefore all non-equilibrium corrections to the distribution functions 
are suppressed. Correspondingly we could start at some initial "moment" Xin ^ 1 
with the equilibrium distribution functions: 

where Tj are the temperatures of the particles, /i is the chemical potential of u^- and 



Ej are the energies, Ej = Jp^ + m? 



At very early times, when ^ f^^vri can neglect the mass of Vr and put 
T^^ = = Ty^ = 1/a = and /i = 0. However, in the case of a large v^- mass 
(e.g. m^^ = 20MeV) we should start at least at Xm = 0.01, in order to satisfy the 
condition ^ m^^. Since the collision rates are very high at this time, the time 
steps in the numerical integration become very small and the program requires much 
CPU time.| Another problem is that a big number of small time steps leads to an 
accumulation of numerical errors. 

Thus it is very difficult to start at ^ m^^ , and we need to find reasonable and 
convenient initial values of Tj and fi. We assume that at > max{10MeV,m,y^) the 
collision rates are high enough to keep all the neutrinos in kinetic equilibrium with 
the distributions given by eqs. (p!^ . 

A simple self-consistent way to find the evolution of Tj and /i in the "almost equi- 
librium" region is to integrate all the kinetic equations over external momenta and to 
get the system of differential equations for the number densities n^^. However, we can 
express the collision integrals in terms of n,y. only in the limit of Boltzmann statistics, 
which differs from the exact Fermi-Dirac one by approximately 10% Therefore 
in order to keep a high precision in the calculations we need another approximation 
for the initial conditions. We require that the temperatures of all neutrino species 
and the photon temperature are equal, T^- = T^, and assume that the initial value 
of the chemical potential ^.u^ is zero. The unknown value of the temperature can be 
found from the energy conservation law (i.e. we do not require that = 1/a). Thus 
our initial conditions at the following: 

f''-^"'^^ = ev/T^ + 1 ' •^^^ = e^^JT^ + 1 ' = e^e/T^ + i ' ^^^^ 

where T^{xin, m^^) is found from the equation p = —3H{p + P), evolved from higher 

temperatures (or smaller x), where expressions (p!5| ) are very accurate with "preinitial" 

'"^CPU time, which is proportional to the number of time steps, grows at least as fast as the rates, 
i.e. as for example, a run with Xin — 0.04 requires twice as much CPU time as a run with 

Xin = 0.05. 
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= 1/a. Note, that for m^^ > IMeV the photon temperature T^{xin,mi,^) strongly 
depends on niu^. 

In order to check the vahdity of our approximation we compared the energy densi- 
ties and number densities for all neutrino species using distributions given by eqs. (|15D 
with the same quantities found from the solution of the kinetic equations (0) for the 
case m,y^ = 20MeV in the time interval 0.04 < x < 0.05. We found that both the 
energy densities and number densities differ in these two approaches only by 0.1%. 

For such values of neutrino mass where m^,^ > Ty( i, non- 

equilibrium corrections cannot be neglected, so for m^,^ > QMeV we choose Xj„ = 0.05, 
and for smaller masses, m^^ < 6MeV, we take Xin = 0.1. 

A careful choice of initial values of Xm and T^{xin-,my^) is very important for the 
numerical results. For example, if for m^^ = 20MeV we take Xin = 0.1 and choose 
to use the "wrong" temperature evolution Ty = 1/a, then we find that the function 
rm = m^^ ny^jn^q is smaller by about 6% than rm calculated with Xj„ = 0.05. With 
our procedure using 7^ 1/a, we are certain that the final results, in particular the 
abundances of light elements, have reached a plateau as long as we use Xi„ < 0.1. 

3.2 Numerical integration of kinetic equations 

For the dimensionless momentum y we took 100 point grid equally spaced in the 
region 0<?/<20. At?/ = the analytical expressions for the collision integrals 
differ from those at y 7^ due to 0/0- uncertainty arising from the factor 1/Eipi in 
front of the integrals and similar vanishing factors in D-functions. Because of that we 
include the point y = separately. This permits to compare the collision integrals 
at y = with those in nearby points, ?/ -C 1, in order to check that numerical errors 
are small in the region of small momentum, y < 1. This is especially important for a 
massive tau-neutrino, because its distribution function rapidly changes in the course 
of evolution. 

Evolution in time is calculated by the simple Euler method for small x {x < 1) 
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and by the Bulirsch-Stoer method (see e.g. ref. IT^) for large x {x > 1). We checked 
that numerical errors in the distribution functions calculated by the Euler method for 
X < 1 are smaller than 0.1%. The Euler method allows us to save a factor of 2 in CPU 
time compared to the second order Runge-Kutta method and more than a factor of 
10 compared to the Bulirsch-Stoer method. For x > 1 the situation is opposite. The 
Euler method requires very small time steps in order to be precise enough, while the 
more powerful Bulirsch-Stoer method requires fewer time steps in order to reach the 
same precision. 

An important check of our program is a run with a small value of z/^ mass {m^^ = 
O.lMeV^) which is in a very close agreement with the results of our calculations in 
the massless case (see ref. [|^). 

4 Discussion of results 

The results of our calculations are presented in several figures below. The energy 
density of massive z/^ is characterized by the quantity rm^^, where r = ny^/riy^ 
and n,^g is the number density of massless equilibrium neutrinos. This quantity is 
presented in fig. 1 for asymptotically large x. Our results are rather close to the 
previous calculations of refs. p, ^, |ll| for small masses, while for large masses they 
are considerably larger than those obtained in refs. O, M and in fine agreement with 



the calculations of ref. . The energy density of z/^ influences the rate of the universe 
cooling and affects two quantities which are important for nucleosynthesis: the frozen 
neutron-to-proton ratio and the time moment corresponding to T ~ 0.065 MeV, 
when formation of light elements begins. These two quantities determine the neutron 
number density at the onset of nucleosynthesis and through that the abundances of 
light elements and in particular of '^He. A larger r results in a larger mass fraction of 
^He both because of a faster cooling rate and because of a more efficient production 
of non-equilibrium electronic neutrinos [^]. Of course the relative abundance of is 
a function of "time" x, so in our numerical calculations we used the exact function 
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r{x) and not just its asymptotic value. The same is true for the energy density, for 
which the exact expression (|^) was used. 

In fig. 2 the product T^a is presented as a function of x for different values of rriy^. 
In the standard calculations this quantity coincides with the ratio T^/T^, which tends 
to 1.40102 for high x. For massless Ur the product T^a is slightly smaller than this 
result because of the energy transfer from the electro-magnetic components of the 
primeval plasma to the massless neutrinos. For a sufficiently heavy Ur the product 
T^a is larger than the canonical value because now the effect of photon heating by 
annihilation of massive Ut is stronger than the energy loss mentioned above. In this 
case massless neutrinos are heated more than electrons, because electrons have to 
share the given energy with photons. For very heavy i/^'s the final photon temperature 
(Ty a) can be estimated from entropy conservation to be (43/36)^/^ times bigger than 
in the massless case. 

In fig. 3 we plot the relative deviations of the energy densities of Ue and from 
the equilibrium value, Sp/p = {p — peq)/peq, where Peq{x) = (7/1207r^)(lMeV"/x)^. 
The sohd and dashed lines in fig. 3 correspond to z/g and respectively. In fig. 3a the 
evolution of 6p/p as functions of x for several different masses of Ur is presented. In 
fig. 3b the asymptotic values of Sp/p as functions of m,^^ are shown. An explanation 
of the different and non-monotonic relative behavior of Sp^^/p and Sp^^/p is given 
below. 

It is convenient to characterize a distortion of neutrino spectrum by the average 
neutrino energy, {E) = p/n. In the equilibrium case it is equal to {Eeq) = 3.15T. In 
fig. 4 we present the relative deviation {{E) — {Eeq))/ (Eeq) for electronic and muonic 
neutrinos. For small masses of Ur this quantity is practically mass independent. In 
this mass range z/g and z/^ are heated by e~e"'"-annihilation. Electronic neutrinos are 
heated more efficiently because of a larger coupling to electrons/positrons due to 
charged current interactions. For higher m^^ the situation is opposite. As we have 
already mentioned, with a heavy Ur, massless neutrinos {ue and z/^) are "overheated". 
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so that their excessive energy is transferred to the electro-magnetic component of the 
primeval plasma (e~ ,7). Now a stronger coupling of v^. to electrons results in a 
more efficient cooling of z/g than z/^. This explains why muonic neutrinos are heated 
less than i/g for a low m^,^, while for a higher vtiy^ the situation is reversed. However, 
for a very large vny^ we return to hotter v^.. This is because very heavy tau neutrinos 
effectively stopped annihilating while i^e, v^l and were still in equilibrium. At that 
stage their temperatures were the same and energy was not transferred either way. 
At smaller temperatures when v^. and z/^ were already almost decoupled, the "old" 
process of e^e^ annihilation heated up z^e more than v^. 

The shape of the electronic neutrino spectrum, disturbed by annihilation, is 
presented in figs. 5 and 6. For large masses the relative distortion 5 is much larger 
than 1 for high momentum, so we present the distorted function itself and compare it 
with the equilibrium one (fig. 5). The magnitude and in particular the non-monotonic 
behavior of the spectral distortions as functions of x for different masses (m = 7 MeV 
and m = 20 MeV) are in agreement with the semi-analytical calculations of refs. 



10| . The effect of the changed neutrino distribution on light element abundance is 
essentially which is plotted as a function of y for different neutrino masses in 

fig. 6. 

In order to extract implications of a massive on light element abundance, we 
have modified the standard nucleosynthesis code (ref. |T5|) in the following way. We 
calculate various quantities at the correct photon temperature and import the values 
to the code. These imported quantities are the neutrino energy densities p^.^ the 6 
weak interaction rates for in ^ ]9)-reactions (this includes reactions (^ and decay 
and inverse decay), and finally (i(ln ar')/dT^ (with the account of neutrinos) which 
governs the evolution of the photon temperature. 

The effect of a nonzero m^^ on the abundances of light elements is illustrated in 
figs. 7 and 8. Comparing our results with the standard calculations of e.g. helium 
production with a variable number of massless neutrino species, we express the impact 
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of massive Vr on nucleosynthesis in terms of the number of equivalent neutrino species 
AA^ = Nf,ff — 3 [^]. The excessive number of neutrino species found from ^ife is 
presented in fig. 7, where our results are compared to those of the previous papers |Q 
|lT| , For low masses we are in a good agreement with ref. 0], while for high 
masses our results are noticeably lower. This difference is probably explained by a 
lower momentum cut-off made in ref. [|lT] . The authors calculated spectral functions of 
different neutrinos up to ymax = '^'^T^/Tu, while we were able to do that up to i/max = 



20. When we diminished our cut-off down to the value of ref. we reproduced 
their higher results in the high mass region. A large sensitivity of ^i/e-production to 
the value of Umax is connected to a strong (quadratic) energy dependence of the cross- 
sections of n ^ p reactions. To independently test the effect of the distortion of the 
high energy neutrinos we calculated the neutron to proton ratio as a function of the 
cut-off of the distortion, numerically solving the kinetic equation for the n/p-ratio. 
The results is, that the n/p ratio is decreasing with an increase of momentum cut-off. 
This continues until rather high momenta (around y = 16 for m^^ = 20) when the 
effect changes sign, so that n/p is increasing when we increase the cut-off in y from 
16 to 20 (as predicted by ref. and refined by ref. ^^). The net result is, that the 
increase of the cut-off from y = 12 to 20 leads to a decrease of n/p in agreement with 
the full nucleosynthesis calculation. With our cut-off in momentum, ymax = 20, the 
error in AA'^i, is approximately an order of magnitude smaller than that at ymax = 12, 
i.e. S (AN^) < 0.02. 

In fig. 8 the number of effective extra neutrino species obtained from other light 
elements: ^H, ^He and '^Li is presented, using 7710 = 3.0. These curves are given for 
illustration only because the best bound on is obtained from '^He. The results are 
in a qualitative agreement with the simplified calculations of ref. |^ , though there are 
quantitative differences. 

The recent confusion with the data on abundance of primordial deuterium (ref. [|1^) 
makes it difficult to put a stringent bound on AA^. It seems that it is rather safe to 
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conclude that AA^ < 1. In this case the consideration of primordial nucleosynthesis 
excludes the mass of Vr in the interval 1 — 22 MeV. Recall that it is valid for the 
sufficiently long-lived Vr-, i-e. for r^^ > 200 sec. Together with the direct experimen- 
tal bound (|l]) it gives m^^ < 1 MeV. This result is obtained for 7710 = 3.0. At lower 
?7io = 1 — 2 the lower bound is slightly strengthend. Hopefully a resolution of the 
observational controversies in the light element abundances will permit to shift this 
limit to even smaller values of m^^. The result in fig. 7 is exact in the sense, that it 
is independent of any measured abundance. One can therefore choose ones favorite 
experiment, extract a c/A^-limit and apply it to fig. 7. In particular, if the limit on 
AA^,^ would return to the "good old" value, AA",^ < 0.3, one could conclude from 
fig. 7 that nii,^ < 0.35 MeV. If one believes in AA^j, < 0.2 (see e.g. ref. |l^) one finds 
m,y^ < 0.25 MeV, and in the case of an over-optimistic limit, AN^, < 0.1, which was 
advocated in the literature a few years ago, we obtain m^^ < 0.15 MeV. 
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A Kinetic equations 

After reducing the collision integral for two-body reactions in eq. (|^) from nine to 
two dimensions the 3 kinetic equations for the distribution functions are solved nu- 
merically. For the sake of brevity we introduce some notation: fa{Pj) = fa \di = 
Di,d2{3,A) = 1)2(3, 4)7^3^4 and 4 = D3/E1E2E3E4, with D-functions defined in 
Appendix A of ref. [|12| . The functions di and are symmetric in all four arguments. 

16 



We introduce one more condensed notation: 



<-(3) f(4). 



(2)^ 



P(l) f(2)/ 



P(3)> 



f(4)> 



and write the coupled system of kinetic equations as follows: 

j dp2P2dp3P3dp4P4S{Ei + E2- E3- E4) 

[6di + 2^2(1, 2) + 2^2(3, 4) - 2^2(1, 4) - 2^2(2, 3) + 6(^3] 



I F(fW f(2) f{3) f{i)\,p(f{l) f{2) r(3) .(4)n 

[di-d2(l,4)-ci2(2,3) + d3] 

I rpffa) f(2) /■(3) /■(4)n , ^(2) .(3) .(4)n 

[2di + ^2(1, 2) + ^2(3, 4) - ^2(1, 4) - ^2(2, 3) + 2^3] 

/i?, fl% & - ^2(1, 3))/E2i?4 

/i?, /i'^) [4(^1 + (c^i - 4(1, 4) - 4(2, 3) + 4) 

+45i5W (4 + 4(1, 2)) /£;3£;4] 

+i^(/i:\ fl'}Ji'^) [^{gl + ^7^) (24 + 4(1, 2) + 4(3, 4) 
-4(1, 4) - 4(2, 3) + 4) -SgL9Rml (4 - 4(1, 3)) /E2E4] } . 

The kinetic equation for f^^ has the same form with the substitutions f^, 

9l- 1- 
The equation for f^^ reads: 

Hxd^cfil^ = J dp2P2dp3P3dpm3{Ei + E2-E3- E4) 

[64 + 24(1, 2) + 24(3, 4) - 24(1, 4) - 24(2, 3) + 64 
+6dimtjEiE2E3E4 + (2(4 - 4(1, 4))/£;2£;3 - (4 + 4(1, '2))/E3E4) 

, rp/^(l) f{2) .(3) , p/ /-(I) f(2) f{3) .(4)y 



/t/„ and 
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[di-d2(l,4)-d2(2,3) + d3- 

-mll2 (c^i + rf2(3,4)/^i^2) 

(di - ^2(1, 4) - ^2(2, 3) + da - yrilj2{d^ + ^2(3, ^))lErE^ 
+4gL9Rml (o?i + 0^2(1, 2) - 2ml^di/ E1E2) / E^E^ 

(2di + d2{l, 2) + o?2(3, 4) - o?2(l, 4) - (i2(2, 3) + 2^3 + m^^(di - ^2(2, 4))/^i^3) 

-S^Lto^e ((c^i - c?2(l, 3)) + 2mldjE,Es) / E^E^ } . 

In these equations the energy 5-function can be trivially integrated away and we 
are left with two-dimensional integral over the energies of incoming particles: 

Hxdjl";^ = / / A{F{f), dk)e{Es + E,-E,- m2)E2pdpmdpA, (17) 



where E2 = E[ + E'^ — Ei, p2 = y E2 — m^, and A{F{f), d^) is the integrand from 
the equations above. 0{E^ + E4 — Ei — 771.2) arises from the energy conservation law. 
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Figure Captions: 

Fig. 1 Relative energy density of a massive tau neutrinos, rm = m^^nu^/riy^, for 
asymptotically large x as a function of m^^ . The solid line corresponds to our results, 
while dashed, dashed-dotted and dotted lines correspond respectively to the results 
ofrefs. iiO]. 

Fig. 2 The quantity T^a as a function of the dimensionless time x for several values 
of the tau neutrino mass in the interval < m^^ < 20MeV. The lowest and upper- 
most solid curves correspond respectively to m^^ = and to m^^ = 20MeV, while 
the dashed curves correspond to several intermediate masses: 4, 6, 8, 10, 12, 15MeV. 

Fig. 3 Relative deviation of the energy densities of massless neutrinos from the 
equilibrium value, Sp/p = {p — Peg)/ Peg, where Peg{x) = (7/1207r^)(lMeV"/x)'^. The 
solid and dashed lines correspond to and z/^ respectively. In fig. 3a the evolution 
of Sp/p with time x for several different masses of u^- is presented. In fig. 3b the 
asymptotic value of Sp/p as a function of m^^ is presented. 

Fig. 4 The relative deviation of the average energy {{E) — {E^q)) / {E^q) for 
electronic and muonic neutrinos as a function of the mass of tau neutrino. The solid 
line corresponds to electron neutrino, while the dashed one corresponds to muon 
neutrino. 

Fig. 5 Spectral distribution of electron neutrinos as function of the dimensionless 
momentum y. Dotted, dashed, and solid lines correspond respectively to equilibrium 
distributions (initial condition), rriy^ = 20MeV and m^^ = 7MeV at asymptotically 
large time. 
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Fig. 6 The change in the electron neutrino spectral distribution times y'^ for 
various masses as a function of momentum y. This quantity is related to the change 
in the neutron to proton ratio, and therefore to the light element abundances. 



Fig. 7 Relative number of equivalent massless neutrino species AiV = N^jj — 3 
as a function of Vr mass, found from ^ife. Our result is presented by the solid line. 
Dashed, dashed-dotted, and dotted lines correspond respectively to the results of 
refs. [ITTl, 



Fig. 8 The number of equivalent massless neutrino species AA^ = iVg// — 3 as 
functions of the tau-neutrino mass, calculated from abundances of deuterium (solid), 
^Li (long dashed), ^He (dashed) and ^He (dotted). 
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Process 


S 




Ve + l^e + 


1/4 


2 [{Pl ■ Pa){P2 ■ Ps) + {Pl ■ P3){P2 ■ Pi) 
+ {Pl ■ P2){P3 ■ Pi)] 


Ve + ^e + Z^iu 


1/4 


\ [(Pl • Pi) {P2 ■ P3) + {Pl ■ Ps) {P2 ■ Pi)] 


Ue + J^e J^T + 


1/4 


\ [(Pl • P4) (P2 • P3) + (Pl • Ps) (P2 • P4) 

- (pi • P2) 


J^e + J^IJ. J^e + ^ti 


1/2 


(Pi • P2) (P3 • P4) + (Pi • P4) (P2 • Pa) 


Ve + J^T J^e + 


1/2 


(Pl • P2)(P3 • P4) + (Pl • P4)(P2 • P3) 
+^^,(Pl •P3) 


Ue + fe e"*" + e~ 


1/2 


2(5i + 5|){(Pi ■P4)(P2 - Pa) 
+(Pi • P3)(P2 • P4)} + AgLgRmlipi ■ P2) 


i/g + i^e + 


1/2 


^{gl + 9r) {{Pl ■ P2){P3 ■ Pi) 

+ (Pi •P4)(P2 - Pa)} - 4:gLgRml{pi ■ p^) 



Table 1: Matrix elements for various electron neutrino processes; Ql = \ + 
siv? 9]v and gji = sin^ 9w Matrix elements for muon neutrino processes are obtained 
by the substitutions Ve — J^i^ and gi ^ Ql — Ql — 
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Process 


S 




l/^ -\- — >• I'r + l^T 


1/4 


2 [{pi ■ Pa) {P2 ■ Ps) + {pi ■ Ps) {P2 ■ Pa) 
+{Pi ■P2){p3-PA) + ^ml^ 

+2m^^ {(Pi • Ps) + (Pi • Pa) - (pi • P2) | 


11^ -\- l/^ — > i/g + i/g 


1/4 


i [(Pi • Pa) (P2 • Ps) + (Pi • Ps) (P2 • Pa) 
- ^^,(P3 •P4)l 




1/4 


\ [(Pl • P4) (P2 ■ Ps) + (Pl ■ Ps) (P2 • P4) 
- "^^,(P3 ■P4)l 




1/2 


(Pl • P2)(PS ■ P4) + (Pl • P4)(P2 • Ps) 

+ml^{p2 -Pa) 


l^T + ^H l^T + l^u. 


1/2 


(Pl • P2)(PS ■ P4) + (Pl ■ P4)(P2 • Ps) 

+m^^(p2 •P4) 


Ur + i^T e"^ + e" 


1/2 


2(^L + 9e) {(Pl • P4)(P2 • Ps) 

+(Pi • Ps) (P2 • P4) - rn\ {p3 ■ Pi) } 
+4^gL9Rml {(pi • P2) - 2ml^} 


Pt + —>■ Vt + 


1/2 


2(^i + ^|){(Pl ■P2)(P3 ■P4) 
+ (Pl ■ P4)(P2 ■ Ps) + "^^,(P2 ■ P4)} 

-^LgRml {(pi ■ P3) + 2m2j 



Table 2: Matrix elements for various tau-neutrino processes; Ql — 91 — ^ — 
— I + sin^ 9w and qr — sin^ 9w 



30 



